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Abstract 

We generalize the free Fermi-gas formulation of certain 3d A/" = 3 supersymmetric 
Chern-Simons-matter theories by allowing Fayet-Iliopoulos couplings as well as mass 
terms for bifundamental matter helds. The resulting partition functions are given by 
simple modihcations of the argument of the Airy function found previously. With 
these extra parameters it is easy to see that mirror-symmetry corresponds to linear 
canonical transformations on the phase space (or operator algebra) of the 1-dimensional 
fermions. 
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1 Introduction and summary 


Supersymmetric localization leads to a dramatic simplification of the calculation of sphere partition 
functions (and some other observables) by reducing the infinite dimensional path integral to a finite 
dimensional matrix model This matrix model can then be solved (sometimes) by a variety of 
old and new techniques to yield exact results. A particular application of this is to check dualities 
— two theories which are equivalent (or flow to the same IR fixed point) should have the same 
partition function. 

In practice, it is often very hard to solve the matrix models exactly, so dualities are checked 
by comparing the matrix models of the two theories and using integral identities to relate them. 
The first beautiful realization of this is in the Alday-Gaiotto-Tachikawa (AGT) correspondence, 
where the matrix models evaluating the partition functions of 4d A/" = 2 theories were shown to 
be essentially identical to correlation functions of Liouville theory as expressed via the conformal 
bootstrap in a specific channel. S'-duality in 4d was then related to the associativity of the OPE 
in Liouville, which is manifested by complicated integral identities for the fusion and braiding 
matrices [3]-[^. 

Here we study 3d supersymmetric theories, which have several types of dualities, of which we 
will consider mirror symmetry and its SLi^^T,) extension [9-14 . Indeed one may use integral 
identities (in the simplest case just the Fourier transform of the sech function) [^ to show that 
the matrix models for certain mirror pairs are equivalent. But is there a way to simplify the 
calculation such that we can rely on a known duality of a model equivalent to the matrix model 
to get the answer without any work, as in the case of AGT? 

Indeed for necklace quiver theories with at least Af = 3 supersymmetry (and one copy of 
each bifundamental field) there is a simple realization of the matrix model in terms of a gas of 
non-interacting fermions in Id with a complicated Hamiltonian 16 . The purpose of this note 


is to point out that the Hamiltonians of pairs of TV = 4 mirror theories are related by a linear 
canonical transformationE] Furthermore we show that the transformations between three known 
mirror theories close to SL{2, Z), which is natural to identify with the S'-duality group of type HB, 
where the three theories have Hanany-Witten brane realizations^] 

In order to demonstrate this we generalize the Fermi-gas formalism of Marino and Putrov to 
theories with nonzero Fayet-Iliopoulos (FI) parameters as well as mass terms for the bi-fundamental 
fields. This is presented in Section where we focus for simplicity on a two-node circular quiver. 

In section [3] we then present the action of mirror symmetry on the density operator of the 
Fermi-gas (the exponential of the Hamiltonian). We also outline the generalization to arbitrary 
circular quivers. The generalization of this formalism to H-quivers and theories with symplectic 


gauge group will be presented in 18 


In the appendix we proceed to evaluate the partition function of the two-node quiver (and its 
mirrors). This was done for the theory without FI terms and bifundamental masses in [^, and 
we here verify that the calculation can be carried through also with these parameters turned on. 
The resulting expressions are not modified much and one can still express them in terms of an 
Airy function. 


16 


^In the specific case of ABJM theory, this was in fact already noted in 

^We should mention of course also the 3d-3d relation 17 , which is closer in spirit to AGT and realizes mirror 


symmetry by geometrical surgery. 


but here we prove it more generally. 
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2 Fermi-gas formalism with masses and Fl-terms 


In this section we review the Fermi-gas formulation 16 of the matrix model of 3d supersymmetric 
held theories and generalize it to a particular A/" = 4 theory that includes all of the ingredients we 
will require for our study of mirror symmetry in the following section. This is a two node quiver 
guage theory with gauge group U{N) x U{N). Each node has a Chern Simons (CS) term with levels 
k and —k. There is a single matter hypermultiplet transforming in the fundamental representation 
of each U{N) factor, and two matter hypermultiplets transforming in the bifundamental and anti¬ 
bifundamental representations of U{N) x U{N). The bifundamental helds have masses mi and m 2 
and each node has a Fayet-Iliopoulos term with parameters and ^ 2 - 

The matrix model for this theory is computed via localisation [^. The result can be easily 
derived by applying the rules presented for instance in pAlIT^I^ 


Z{N) = 


{N\y 




ni<,- 4sinh^ 7r(Af ^ - X)'’)4:smh^ 7i{Xy’ - X)-"’) 


Kj 


(l)^ 

J 


.(2) 


,(2)^ 


Hij 2 cosh 7r(Af ^ - Xf'’ + mi) 2 cosh 7r(Af^ m 2 ) 


^ g27riCiAT+7rife(Al^*)2g27riC2A)^^-7rjfc(AP^)^ 

Xll- 

2=1 


( 2 . 1 ) 


2 cosh ttA,-^^ 2 cosh ttA,-^^ 


The crucial step in rewriting this expression as a Fermi-gas partition function is the use of the 
Cauchy determinant identity 


- Vj) 

UiA^i - Vj) 


N 




ctGSn 


2=1 


2/cr(*)) 


( 2 . 2 ) 


Applying this to (2.1) we may write the partition function as 

N 


Z{N) 


(Niy 


ctiS 5 jv *=1 


*=1 2cosh7r(A) 


( 1 ) 




N 


"" SJ" E[2cosh7r(Af^-A^'^).)+m2)n 


g27riCiAh^-|-7rifc(A)^')2g27riC2Ap^-7rjA:(Al^^)2 


2 cosh :rAy 2 cosh n'Ap* 


(2.3) 

( 2 ) ( 2 ) 

A relabelling of eigenvalues A) —)• allows us to resolve one of the sums over permutations, 

pulling out an overall factor of A^! giving 


= m 


5^(-l)” f ci«A<‘'d«Af> n 

crG5'iv 


jA g27r<iA)^^gj7rfc(AA)2 


2 cosh ttA-^^ 2 cosh7r(A— A^^ -|- mi) 

1 


g 27 riC 2 Af'g-i 7 rfe(Af 


1 

TV! 


2 cosh vrAf ^ 2 cosh 7r(Af ^ -|- m 2 ) 

f5:(-l)'/d-A<‘'/A(A<‘>.A“). 


(2.4) 


ctSSjv 
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Here we expressed the interaction between the eigenvalnes in terms of the kernel K, which 
can be considered the matrix element of the density operator K dehned by 

p2TTiCliq-\-'Kikq^ p2nimip p27rim2P 

K{quq2) = {qi\K\q2), K= —— . . (2-5) 

2 cosh Tiq 2 cosh Tip 2 cosh irq 2 cosh Tip 

where p and q are canonical conjngate variables [q,p] = ih with h = l/27r and we have made nse 
of the elementary identities 


f{q) \q) = f{q) \q) 

^-2nimpj^^ynimp ^ j 

tel—h^k2>= ' 


cosh Tip 


cosh7r(gi - ^ 2 ) 


( 2 . 6 ) 

( 2 . 7 ) 

( 2 . 8 ) 


To stndy the system in a semiclassical expansion it is nsefnl to represent the operators in 
Wigner’s phase space, where the Wigner transform of an operator A is defined as 


Aw{q,p) = I dq' {q-^ 


A 




Some important properties are 

(yAB)\\r = Ayi/ ic By/, ic = exp 


ih 

~2 




Tr(i) = 


dqdp 

2Tih 


A 


w ■ 


(2.9) 


( 2 . 10 ) 


For a detailed discnssion of the phase space approach to Fermi-gasses see 16 and the original 


22 


paper 21 . For a more general review of Wigner’s phase space and also many original papers 
see 


In the langnage of phase space the kernel K (2.5) becomes 

p27ri<^iq-l-7rikq^ p27rimip p27Ti(^2q—7T'i’kq^ 

Kw = 


^27Tim2p 


2 cosh Tiq 2 cosh Tip 2 cosh Tiq 2 cosh Tip 


( 2 . 11 ) 


Clearly the partition fnnction can be determined from the spectrnm of K or Ky/- The leading 
classical part comes from replacing the star prodnct with a regnlar prodnct. In the appendix we 


ontline the calculation of the partition function, extending 16 


3 Mirror symmetry 


In this section we examine the theory studied in the previous section, with vanishing CS levels 
(see top quiver in Figure [^, where the density function (2.11) becomes 


:,2mCiq 


^2'Kimip 


o2.mC.2q 


^2nim2P 


Kw = 


* 


* 


(3.1) 


2 cosh Tiq 2 cosh Tip 2 cosh Tiq 2 cosh Tip 
It has been known for a long time that this theory has two mirror theories, related in the IIB brane 


construction by S'L(2,Z) transformations 12 . As we show, the density functions of these theories 
are simply related by linear canonical transformations. 
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Figure 1: Quiver diagrams summarising the two node theory we discuss in the text and its two mirror 
duals. Each circle represents a U{N) vector multiplet, labelled inside by the CS level k and ontside 
by the FI parameter. Edges represent hypermultiplets. Those connecting two circles are bifundamental 
fields with the mass indicated next to them. The boxes represent U{1) flavor symmetries of fundamental 
hypers, which in our examples are massless. 


3.1 S transformation 

The first of the known mirror theories is one with identical matter content but with mass and FI 


parameters exchanged 11 


mi -)■ mi = -Cl, m 2 -)■ m 2 = -C 2 , Ci -t Ci = "^ 2 , C 2 -t C 2 = • (3.2) 

This is illustrated by the bottom right quiver in Figure 
At the level of the density function, this gives 


g—27ri^ip ^27vim\q ^—2'Ki(^2P ^27vira2q 

-k -:- -k -:- k: 


^27Tim2q g-27r<ip ^27rimiq g-27r<2P 

^ ^ ^ _ 

^ 2 cosh ng 2 cosh np 2 cosh vrg 2 cosh np 2 cosh irp 2 cosh nq 2 cosh np 2 cosh Trg ’ 

(3.3) 

where the last relation represents equivalence under conjugating by 2 coshnq ' 
density is the same as (3.1) under the replacement 

p ^ q , q ^ —p . 


(3.4) 


3.2 U transformation 


To get the second mirror theory we apply to (3.1) the replacement 


p ^ p + q , g —)■ —p. 


(3.5) 
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The result iE 


= 

^ 2cosh7rp 2cosh7r(p + g) 2cosh7rp 2cosh7r(p + g) 


g-27r<ip ^2nim-i{p+q) 


g-27r<2P 
-k -:- -k 


^2mm2(p+q) 


g-27riCip 

2 cosh 7Tp 


k e 


—iTrq 


^2'Kimip 

2 cosh Tip 


k e 


iTrq 


^-2wiC2P 
2 cosh Tip 


k e 


—iTTq 


^27vim2P 

2 cosh Tip 


(3.6) 


^ e 


tnq 


In the second line we have made use of the identity 


e-"’ *f(p)*e-'> =f(p + q). 


(3.7) 


To read off the corresponding quiver theory from (3.6), each +2Cg) term can be associated to 
a U(N) node with CS level k and FI parameter C,, while each 2 coshnp comes from a bifundamental 
hypermultiplet with mass m. The transformed density operator corresponds therefore to a circular 
quiver with four nodes that have alternating Chern-Simons levels k = ±1 and vanishing FI param¬ 
eters. The bifundamental multiplets connecting adjacent nodes have masses {—(i, mi, —( 2 ,^ 2 }, 
as in the bottom left diagram of Figure [T] 


As we discuss in the appendix, the partition function can be expressed in terms of Tr (A.l), 


which in phase space is given by an integral over p,q (2.10). Since we showed that mirror symmetry 


can be viewed just as a linear canonical transformation, which is a change of variables with unit 
Jacobian, it is clear that mirror symmetry preserves the partition function. 

3.3 SL{2,Z) 

It is easy to see that the transformations we used in the previous sections close onto SL{2,Z). 
Indeed, dehning T = SU we find the dehning relations 


S^ = -I, {STf = I. 

More general SL{2, Z) transformations will give density operators with terms of the form 


(3.8) 


cosh7r(ap -|- bq) 


(3.9) 


The cases with a = 0 and 6=lora = ±l and 6 G Z have a natural interpretation as a contribution 


of a fundamental field, or as we have seen in (3.7), from conjugating the usual 1/coshvrp by 
CS terms. But these manipulations cannot undo expressions one finds from a general S'L(2,Z) 
transformation of Kw- In these more general cases, the transformed density operator can still be 
associated to a matrix model, but it cannot be derived from any known 3d lagrangian. 

This is also manifested in the IIB brane realization, where any SL{2, Z) transformation will 
lead to some conhguration of {p, q) branes. Most of those do not have a known Lagrangian 
but one could associate to them a matrix model 


description 23 


14 , which would indeed lead to 


the transformed density operator. 


particular under (3.5). 


^Note that the definition of the star product (2.10) is invariant under linear canonical transformations, and in 
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3.4 Mirror symmetry for generic circular quiver 

The manifestation of mirror symmetry as a canonical transformation naturally generalises to the 
entire family of A/” = 4 circular quivers with an arbitrary number of nodes. Applying the Fermi-gas 
formalism, it is easy to see that the density function for such a theory with n nodes is given bj|^ 


Kw — 


^2TTi(aq 


^2mmaV 


* (2 cosh Tiq )2 cosh up ’ 


(3.10) 


where C,a denotes the FI parameter of the node, Na denotes the number of fundamental matter 
fields attached to the node and rria denotes the mass of the bifundamental field connecting the 
and (a + 1)*^ nodes. 

We can now apply the S and U transformations of the previous section, and look to see if 
the resulting density functions can again be interpreted as coming from the mirror gauge theories. 
Applying the S transformation we get 


K 


(s) 

w 


n 


g-27riCaP 

* (2 cosh 7rp)'^“ 


^2mmaq 

2 cosh nq 


(3.11) 


This density function is that of a circular quiver theory with X]a=i nodes and n fundamental 
matter helds. The fundamentals are attached to nodes which have FI parameters rria-, and are 
separated by A^^ — 1 other nodes. The masses of the bifundamentals connecting them add up to 

-C.0 

Applying the U transformation we get 


AT=n. 


^-2wiCaP 

\Na 


^2vima{p+q) 


n 


g-27riCaP 


=1 * (2 cosh vrp) “ 2cosh7r(p + g) * (2 cosh vrp) 


\Na 




^2mmaP 

2 cosh np 


-k e 


mq 


(3.12) 


The mirror theory can be readily read off from this density function as a circular quiver theory with 
X]a=i + n nodes and no fundamental matter. Each node has Chern-Simons level k = +1,-1 
or 0. Further details concerning the mass parameters and value of the Chern-Simons level at each 
node can be read off in much the same way as for the previous example. 

A further generalisation we have not yet considered is to turn on masses for the fundamental 
fields. This corresponds to replacing each of the (2cosh7rg)“^“ in (3.10) with a product of Na 
terms with masses pi 


Na 




^2mCaq ^ 


*^-^*2cosh7r(g -|- pi) 


Na 

= * Yl 

i=l 


^2mpip ^ 


2 cosh 71 q 


* e 


— 2nifiip 


_ g-27r*Ca/ilg27ri/iip ^ 


„27riCag 


2 cosh Tiq 


Na 


^ g-2«/.ip ^ J-J /g2w^ 


i=2 


2 cosh Tiq 


•k e 


—2nipip 


(3.13) 


'^The product is defined by ordered star multiplication. 

® At the level of the matrix model, this additional freedom to choose mass parameters in the mirror theory simply 
amounts to the freedom to make constant shifts in the integration variables. 
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Where in the second line we chose to associate the FI term to the hrst fundamental held, picking 
up an overall phasej^ 


Once we apply S' or f/ transformations to (3.13) it becomes clear that these mass terms become 
additional FI parameters, as is expected. 
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A From density to Airy function 

In this appendix we outline the computation of the large N partition function for a theory with 


mass and FI parameters. Since the calculation follows closely the method outlined in 16 we will 


be rather cursory and refer the reader to 16 for more detail, highlighting the new features that 
appear due to the FI and mass parameters. 


To evaluate (2.4), one notices that it combines to give products of Zi = Tt^KwY where the 
values of I depends only on the conjugacy class of a. Instead of summing over all permutations we 
can sum over conjugacy classes which have mi cycles of length I and the combinatorics give 


z(a') = ETI 


{mi} I 


miU^i 


(A.l) 


where the primed sum denotes the restriction to sets that satisfy YYi Following the usual 

analysis from statistical mechanics 24 we consider the grand canonical partition function given 
by 


l(l^) = 1 + 5^ Z{N)e>^^ = Exp 


N=1 


-E 

. i=i 


I 


(A.2) 


We consider the density function (3.1), and using (2.7), (2.10) rewrite it a^ 


Kiy = e^*(”^U2-|Cimi-|C2m2)g7rj((2Ci+C2)<?+mip) ^ 


2cosh7r(g+^) 2cosh7r(p—y) 


2cosh7r(g —^) 2cosh7r(p+y) 

^ ■n-i((2g+{'mi+2m2)p) 


(A.3) 


-k e 


®There is a freedom to distribute the FI terms arbitrarily among the fundamental fields, leading to a different 
phase in front, see also footnote]^ 

^It is also possible to rearrange the expressions such that one p is not shifted and the other shifted by C 2 and 
one q is not shifted and the other shifted by mi. We choose this more symmetric expression for later convenience. 
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Figure 2: Fermi surfaces for F = 20 and masses m = 0 (a) and m = 1 (b). The red lines are the exact 
classical Hamiltonians (A.8) and black lines the polygon approximat ion (|A.9 ). The region between them 
is the leading order perturbative correction given in the first line of (A. 13). 


In order to get a hermetian Hamiltonian below, we specialize to the case 

Cl = -C 2 = C , mi = -m 2 = m , (A.4) 


and conjugate 

Kw (2 cosh 7r(g - f)) U Kw * 

This gives the kernel 

g-27riCm 


(2cosh7r(g — yjja 


(A.5) 


Kw = 


(2cosh7r(q'- y ))2 2cosh7r(p-|) 2 cosh7r(q'+ y) 2cosh7r(p+|) (2cosh7r(g - y ))2 

(A.6) 

The phase in Kw leads to an overall phase (e front of the partition function and can 

be removed and reintroduced at the end ( A.lQp . 

Following Section 4 of we compute the partition function by studying the spectrum of the 
one particle Hamiltonian 

Hw = —^og^Kw- (A.7) 

To hnd an expression for Hw one must perform a Baker-Campbell-Hausdorff (BCH) expansion of 
the logarithm in (A.7). Setting C = Oj^the leading classical term in this expansion is simply 


m 


Hc\ = log 2 cosh TT g + — ) ) + log 2 cosh tt g — — ) ) + 2 log(2 cosh up) 


m 


(A.8) 


The effects of nonzero C and nonzero m are completely analogous. 


























For large p, q this is 


H, 


cl 


TT 


m 

~2 


TT 


m 

~2 


2vr|p| 


(A.9) 


It is clear that the approximate Hamiltonian is independent of m for |g| > y. In Figure we 
display the exact classical Fermi surface and polygonal approximation for a particular value of E 
and with vanishing and non vanishing mass. The only change to the polygon from turning on the 
mass is the removal of the two triangles with bl > ^ ~ y whose combined area is m?/2. 

The number of states below the energy E is given by the area enclosed by the curve H = E. 
Using the polynomial approximation this is just 


n{E) 


dqdp 9[E — TT 


m 

~2 


TT 


m 

~2 


) Tn^ 1 

= CE^--, C = ^. (A.IO) 


This expression is only approximate and gets corrected by accounting for the difference between 
(A.9) and the exact quantum Hamiltonian (A.7). We do this by modifying the number of states 
to 

(A.ll) 


m 


n{E) = CE^ - — + no + n^piE ), 

We outline the calculation of no below. The main point is that it does not depend on m. n^p denote 
nonperturbative, exponentially suppressed corrections at large but are such that n(0) = 0. To 
satisfy this, n^p clearly has to depend on m, but this will have no effect on our end result where 
we ignore the nonperturbative terms. The approximation ( |A.9 ) is valid where both p or g are 
large, and as shown in 16 , the quantum corrections to (A.8[) (from the BCH expansion of (A.7)) 


are also exponentially suppressed there. All the corrections are therefore associated with regions 
where either p or g are small, namely around the vertices of the polygon. We then consider 
the contributions to np from each region separately, integrating in each case the (perturbative) 
corrections to the boundary. For instance, around g ~ 0, p 3> 1 the hrst quantum corrections of 
the Hamiltonian are given (up to terms exponentially suppressed in p) by 


vr 


Hcl —)■ + — 


cl 


24 \ coslF 7r(g + f) 


2 

cosh^ 7r(g 



(A. 12) 


The difference in the area between the polygon and the quantum corrected Fermi surface for large 
E approaches 


r ^ 


m 

+ vr 

m 

1-00 2vr 

vr 

g H- 

^ 2 

CN 

1 


TT^ / 1 

24 (^coshV(g + f) 


log 2 cosh vr ( g + — 


m 


log 2 cosh vr (g 


_2_t 

cosh^ vr(g — y)y 


1 

^ ■ 



(A.13) 


As advertised, this is independent of m (which can be seen by splitting the integral into two term 
with g ± m/2 and shifting the integration variable). 


®These also include the so called Wigner-Kirkwood corrections to the formula (A.IO), which manifest as boundary 


integrals on the Fermi surface and are nonperturbative for the same reason as in 


16 
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The analog expression aronnd the p = 0, g S> 1 vertex is 


dp f TT^ 

/ — 27r|p|-21og2cosh7rp-—-^— 

J_^ 2tt \ 48 cosh Tip 

Snmming over the contribntions from all four regions we get 

^ / 1 5 \ 7 

~ 48 J “ ■ 
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(A. 14) 


(A.15) 


It is not hard to see that all the higher order quantum corrections do not modify riQ and in 


particular do not depend on m. See the discussion in Section 5.3 of 16 


From n(E) it is easy to calculate the matrix model partition function. The grand canonical 

dn{E) 


potential, the logarithm of (A.2) is 


J(/i) = logS(/i) = / dEp{E) log(l + ^), p{E) = 


dE 


(A.16) 


At large p this integral reduces to 

C 


d (h) ~ + h ( ^0 ^—I —^ ) + A + 0{pe . 


A is a constant that we will not concern ourselves with (it is studied in more detail in 
subsequent papers). From this we can extract the canonical partition functioij^ 


Z{N) = — / = C-h^ki 

27ii 


C 3 ( iV — no- - —h — 


+ 0{e 


-N\ 


(A.17) 
and 

(A.18) 


It is straight forward to include also the FI parameter Remembering the extra phase in (A.6) 
one hnds 

7i‘^C^m^ + e 

C ^ { N — no ■ 


Z{N) = ki 


+ 


(A.19) 


3 2 

One can treat quite general kf = 3 necklace quiver theories in a similar fashion, subject to the 
technical constraint that the sum over FI parameters and the sum over bifundamental masses both 
vanish. The analog of (A.9) will again be a piecewise linear Hamiltonian 

H ^'^\aiq + biP +Ci\, (A.20) 


The parameters Oj and 6* are determined by the Chern-Simons terms and q are due to mass and 
FI terms. The volume of the corresponding polygonal Fermi surface is again of the form 


CE^ + B 


(A.21) 


Similar arguments to those made above guarantee that the full c* dependence appears via a shift 
B which can be found already in the polygonal approximation. 
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For a discussion of the integration contour, see 


25 
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